Temperature effects on mixed-state geomètric phase 
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Geomètric phase of an open quantum system that is interacting with a thermal environment (bath) is studied 
through some simple examples. The system is considered to be a simple spin-half particle which is weakly 
coupled to the bath. It is seen that even in this regime the geomètric phase can vary with temperature. In addition, 
we also consider the system under an adiabatically time-varying magnètic field which is weakly coupled to the 
bath. An important feature of this model is that it reveals existence of a temperature-scale in which adiabaticity 
condition is preserved and beyond which the geomètric phase is varying quite rapidly with temperature. This 
temperature is exactly the one in which the geomètric phase vanishes. This analysis has some implications in 
realistic implementations of geomètric quantum computation. 

PACS numbers: 03.65.Vf, 03.65.Yz 



I. INTRODUCTION 

Geomètric phase (GP) of a quantum state is one of impor- 
tant concepts in quantum theory. Indeed, Berry was the first 
who emphasized the geomètric nature of the phase acquired 
by an eigenstate of an adiabatically-varying Hamiltonian in a 
closed loop in parameter space []]]. Since its discovery, this 
concept has been a subject of a vast investigation and gen- 
eralization in many aspects 0. An important reason for 
the interest in the concept of GP is its relevance to geomètric 
quantum computation [4]. Indeed, it is believed that the purely 
geomètric characteristic of such phases potentially provides 
robustness against certain sources of noise 01 . 

For general evolutions of an open quantum system, where 
the dynamics is generally nonunitar y, th e GP has been de- 
fined in different methods 000 0M El Cl 13 E Cl 
The first general approach is Uhlmann's mathematical method 
which is based on a purification of mixed states 0. Another, 
more recent, approach is based on a kinematic extension of 
the interferometric approach firstly used in the case of unitary 
evolutions |jl3ll. These methods are generally argued to be 
different 11611 . however, recently a formal approach bridging 
between them has also been proposed fÏ7ll . 

As mentioned above studying how GP is affected by en- 
vironmental noises is an important issue in investigation of 
robustness of geomètric quantum computation. In some ear- 
lier works, effect of different t ype s of decoherence sources on 
GP has been studied. In Ref. II 211 . it has been shown that the 
GP of a spin-half system driven by one or two mode quantum 
fields subject to decoherence behaves differ ently for adiabatic 
and nonadiabatic evolutions. Also in Ref. 11811 . the GP of a 
two-level system driven by a quantized magnètic field sub- 
ject to dephasing has been considered and it has been shown 
that the GP acquired b y th e system can be observed even in 
long time-scales. Ref. 11411 addresses the change of the mean 
GP when the two-level system is weakly coupled to a thermal. 
An interesting feature of this investigation is that the mean GP 
does not have any thermal correction up to the first order in the 
coupling constant, that is, any dependence on temperature is 
of higher orders in the coupling constant. In Ref. 11911 . it has 
been shown that coupling to an environment induces some ge- 



omètric and nongeometric corrections to the Berry phase (BP) 
of a spin-half system that is under an adiabatically slow rotat - 
ing magnètic field. Moreover, it has been argued that the BP 
can be observed only in experiments whose time-scales, on 
one hand, are slow enough to ignore nonadiabatic correction, 
and fast enough, on the other hand, to prevent dephasing from 
deleting all phase information. Existence of such a finite (adi- 
abaticity) time-scale for the GP, indeed, has been shown to be 
a general feature of open quantum systems fÏ5l Eoll . Some 
other related studies about decoherence effects on the GP can 
be found in Refs. EHHHHÜ. 



II. THE MODELS 

In this paper we consider a simple spin-half particle in a 
magnètic field B(t) which is weakly coupled to an environ- 
ment (a thermal bath, for example, consisting of photons) 
with temperature T. We are mainly interested in finding ther- 
mal effects of the environment on the GP of the system. To 
be specific, here, we follow the kinematic definition of the 
GP JÏ3ll . If the density matrix of the system evolves as 
0(0) -> g(t) = EiPi(t)\wi(t))(wi(t)\, where t e [0,r] and 
Pi(t)'s (\wi(t))'s) are eigenvalues (eigenvectors) of g(t), then 
the GP acquired by this state during this interval reads as fol- 
lows: 

$ = argE i [Pi(0)Pi(r)]5(w i (0)| Wi (r))e-/o T <- í (*)l^(*)>dt. (1) 

Throughout this paper we always consider the GP measured 
in mod 2it. To use the above relation one must solve the 
eigenvalue problem for the density matrix of the system to 
find eigenvalues as well as the corresponding eigenvectors. 

Before discussion about coupling of an environment and 
the system in a relatively realistic model, let us illustrate a 
couple of simple examples that can reveal some general fea- 
tures of more complicated models. As the first example, sup- 
pose that due to interaction with the environment, our sys- 
tem which is a spin-half particle, initially is prepared in the 
thermally impure state (1 — e)|+; 0} n (+; 0| + eg±., where 

/ cos Èlíl e " í v{t) \ 

|+;t)„ = JZ , (t) is an eigenvector of n.(7, n = 



FIG. 1: GP ($, in radian) vs. temperature (T, in the units where 
h = k,B = 1) for different vàlues of Here, we have taken e = 1/3 
and A = 1 (in the same natural units). 

(sin$cos</?, sintfcosy;, costf), g th . = A| + ;0)„(+;0| + (1 — 
A)| — ; 0) n (— ; 0|, y^j = e~ A / T (we work in the natural units 
where h = ks = 1), A is the energy gap of the system, 
a = (a x , a y ,(J z ) is the vector of the Pauli matrices, and for 
simplicity we take $(í) = i?(0) and ip(t) = ujat. Now let us 
assume that the density matrix of the system has the following 
(unitary) evolution during the period t € [0, ^] 

[(1 - 6) + Ae]|+;t)„(+;í| + (1 - A)e|-;í)„(-;í|. (2) 

By using the general formula one can easily obtain the GP 
acquired by this state after its cyclic evolution as follows: 

$ = 7T + arctan(tan(7rcostf)[l - 2e/(l + e~ A / T )]) .(3) 

Figure[2shows that depending on the value of -d the GP may 
decrease or increase with temperature. We will see a bit later 
that the GP shows such a typical behavior even in more com- 
plicated cases. 

As the second example, let us consider two coupled spin- 
half partides, one of which driven by a time-varying magnètic 
field, with the Hamiltonian 

H(t) = \B(t).ax + J(a +1 a_ 2 + H.c), (4) 

where a± = \{u x ± ia y ) and the indices 1 and 2 indicate 
the corresponding subsystems. The magnètic field is consid- 
ered to be Bn(t), where n(í) has a form as in the previous 
example but it rotates about the z-axis adiabatically. In fact, 
in this simple case the second particle can be taken as an en- 
vironment which is coupled to our system of interest. The 
eigenvalues (Ej) and instantaneous eigenvectors (^(í))) of 
this Hamiltonian can be found simply as follows: 

E 1{3) = ~E m = ^l + ^±f N A 2 + 4sin 2 í ?, (5) 
= -^=(^-1 U) + b 3 \ TT) + Cj\ U) + dj\ |T)), (6) 

where d, = "^""t^ e^*, ò 7 = - g ' +c ° s V e-'""*, 
cij = smiïe~ luJ °\ Cj = Ej - cosí?, N 3 = \a 3 \ 2 + b) + 

+ \dj\ 2 , and k = ^. In Ref. |Q, the BP's of the subsys- 
tems for all eigenvectors have been calculated, and it has been 
shown that the BP of the composite system is sum of the BP's 
of the two constituting subsystems. Here, we suppose that 




FIG. 2: BP ($, in radian) of the subsystem 1 vs. temperature (T, in 
the units where h = ks = 1), obtained from p[, j = 1, 2, 3, 4, for 
-ff = 7r/4, with e = 0.1. The numbers show the related vàlues of the 
coupling constant k. 

the composite system is in the thermally diluted state p 3 ~ 

(l-e^jX^I+epthermai, in which Ahemiai oc e~ H/T . We are 
interested in calculating the BP of the subsystem 1 to see how 
it is affected by this simple thermalization mechanism. By a 
simple calculation, it can be seen that in the Schmidt decom- 
position of the eigenvectors coefficients are time-independent, 

= E a \frt\FÍ(t))i\fí(t))2- This fact in turn re- 
sults into the reduced density matrix p{(í) = tizipj) = 
12ia Pa\Fa(t)) (Fa(t) \ for the subsystem 1, in which for i = j 
we have p a = [1 — e + ee~ Ej ^ T /tr(e~ H ^ T )]p 3 a and p l a = 
te~ Ei l T p l a / ü{e~ H l T ) otherwise. To calculate the BP of this 
density matrix we use the idea of purification 

*Ben-y(pj(í)) = $Beny (|T* (í))l a ) , (7) 

where |TJ(í)) la = ^ a ^I^(*))il4>a ^ apurified state 
obtained from p\ (t) by attaching an ancilla a. By using this 
purification and considering constancy of the coefficients of 
the Schmidt decompositions, it can easily be concluded that 

^Beny( / 0Í) = E i aPÍ $ B et ry(^) J (8) 

where $ B erry(^) = i ffi {F^ (í) \í* (í))dí. Figure illus- 
trates the BP of the subsystem 1 during a cyclic evolution, for 
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FIG 3: BP ($, in radian), vs. temperature (T, in the units where 
h = kB = 1) for different vàlues of e. Here, we have taken p\, and 
k = 2, í? = n/4. 
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FIG 4: GP ($, in radian), vs. mean number of photons (n) in the 
environment. Here, wo = 2, $ = ^ (a), f + t (b)> and the num- 
bers associated to each curve denote the related value of the coupling 
constant k. 



j = 1, 2, 3, 4. As is seen, after a decrease for some temper- 
atura ranges, when temperature increases the BP increases as 
well. Figure[3]shows that the more thermal the composite sys- 
tem is the lower the value of the BP is. However, even in this 
case the general behavior of the BP vs. temperature is as in 
Fíg.El 

Now, let us consider a more realistic case in which the spin 
is interacting with a thermal environment consisting of pho- 
tons. It is known that when the spin is not isolated its dy- 
namics is much more involved. When the magnètic field is 
not time-varying, interaction with the environment generally 
induces energy and phase relaxation processes (with the time- 
scales T\ and T 2 , respectively), in addition to a Lamb shift 
of the energy levels (SE^^) 12611 . It is natural to expect that 
when dephasing is the leading regime of the interaction no 
GP effect could be observed |19]. However, in the weak- 
coupling limit, defined by B T 2 , in experiment time- 
scale r these effects can be seen if both adiabaticity condition 
B ^ t~ 1 and no-dephasing regime condition r < T 2 can 
be satisfied simultaneously IÍ19I1 . To fulfill these conditions, in 
our discussion we assume that the field varies adiabatically 
on a closed loop in the parameter space in time r. Addi- 
tionally, we also suppose that the time-scale in which typ- 
ical correlation functions of the environment decay, í env ., is 
much smaller than r and all dissipation-induced time-scales, 
i.e., í env . < r,T<ji ss . = mm{T 1 ,T 2 ,5E~J [Bb }. These assump- 
tions overall enable us to use Markovian approximation in the 
form of a Lindblad-type master equation [27] 

q = -i[H, o\ + \Y.k (IT*, <F\] + [Tfce, r£]) , (9) 

where H is the effective Hamiltonian of the (open) system 
and Tfc's are the quantum jump operators describing system- 
environment coupling. The magnètic field is supposed to be 
in the direction n(í) = (sini? cos ^{t),smd sva.Lp{t), cosí?), 
where <p(t) = wqí, and r = ^j- is the time necessary for 
completion of a loop. The Hamiltonian of the system is 
H(t) = ÍB(t).a, where a z = |e)(e| - \g){g\ = (J^), 
in which \g) (|e)) is the ground (excited) state of the system in 
the absence of external magnètic field. For such an open sys- 
tem the quantum jump operators are Ti = \J n{fi + ï)|<j)(e|, 
which describes spontaneous and stimulated emissions into 
the bath, and T 2 = \/~Kïï\e) {g\, which is the corresponding 
absorption from the bath, where n(u>o, T) = (e"°/ T — is 




FIG. 5: The path traversed by the tip of the Bloch vector of the spin 
starting from & = f , forn = 0, 1, 5, 10, 20, 80 (left), and = vt/4, 
for n = 0,1, 5, 10, 20 (right). In both of the figures, we have set 
k = 10 -2 , B = 10 3 , and lu = 2. 



the mean number of photons in the bath. A master equation in 
the Lindblad form usually can be applied in the case of time- 
independent Hamiltonians under Markovian conditions. The 
simpler case in which the Hamiltonian of the system is ^ujo<j z , 
a constant magnètic field along z-axis, has been studied ear- 
lier in Ref. Í14ll . and to the first order in coupling constant k, 
the mean GP has been calculated. It has been shown that up to 
this order no thermal correction appears in the GP, that is, any 
dependence on temperature (through n) is of higher orders of 
K. Figure|4]depicts temperature-dependence of the GP for the 
two vàlues of d = ir/2, 3ir/4 and different vàlues of the cou- 
pling constant k. For the vàlues i? < f the case is similar to 
$ = \. When •& ^ 5, this figure shows a decrease in the GP 
with temperature, whereas for -Q > f the GP can also show 
an increase for some temperature range whose size decreases 
when k increases. 

Now, let us consider the general case with a time-varying 
magnètic field, as explained above. To be able to use the 
Lindblad master equation we change our reference frame 
and go to the diagonal frame in which the Hamiltonian 
has the simple time-independent form ha z , where h = 

ii?^/sin 2 i?+ (cosí? - u> /B) 2 Í2T1 . Here it is worth not- 
ing that another alternative method to consider the problem 
is to define an adiabatic Hamiltonian I2ÜI . which in the ab- 
sence of interaction with environment shows the BP effect. In 
this way, in contrast to the rotating frame method, only the 
length of magnètic field changes. The results of this method 
have been shown to be equivalent to the rotating frame method 
in the first order of the rotation frequency 12311 . The den- 
sity matrix of the system in the diagonal frame, denoted 
by g, is related to the density matrix in the lab frame as 

$ 7r 
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FIG 6: GP ($, in radian) vs. mean number of photons (n) in the 
environment. Here, wo = 2, k — 10 -3 , B — 10 3 , and the azimuthal 
angles d = f ; ï > ï + 5 are chosen for this plot. 
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FIG. 7: (a) GP ($, in radian) vs. mean number of photons (n) in 
environment (thin curve) and W = |(wi(0)|wi(t))| (thick curve). 
Here, uj a = 2, B = 10 3 , k = 5 x 1(T 3 , = \ . (b) n ad . vs. 
coupling constant k. Here, we have set cjq — 2 and £? = 10 3 . 



g(í) = e-% tWoa *VQ(t)Ve% tUoC 
and a± 



where V = a_a x + a + a z 



| ± ^(-B cosí? — a>o). The explícit form of the 
master equation in this frame is as follows: 



g =—i[ha z , g] + n(n + l)(2er_£Í(7 + — ga + <j- — cr + <j-g 
+ KÏï(2a+g(T- — gtj-(j + — <j-cr + g). 



10) 



For simplicity we have omitted the Lamb shift terms, because 
the Lamb shift contributes to the GP for r ^ 2tt/lüq only Q. 
These equations can be solved exactly as follows: 



GP is very tedious, instead, we plot the results of some special 
cases in Fig. 0. It shows that for a relatively large temper- 
ature range the GP is varying very slowly with temperature. 
This is in accordance with the results of Ref. 1 14|. In addition, 
Fig. 0a) shows that when temperature gradually increases the 
GP decreases so that in some threshold value it vanishes. After 
this value the GP shows a rapidly varying behavior. The value 
of this threshold depends on the parameters of the problem. 
An special and important feature of this figure is that the over- 
lapof |* ) = |wi(0)) and |wi(t)), W = |(wi(0)|wi(r))|, 
as a measure of adiabaticity of the evolution vanishes in the 
same threshold temperature. This indicates that after this tem- 
perature adiabaticity condition breaks down. This implies ex- 
istence of a temperature-scale n^. for the observation of the 
adiabatic GP. This is similar to the existence of a finite adi- 
abaticity time-scale for observation of the GP, as pointed out 
in Ref. ÍI5I1 . In fact, this similarity between temperature and 
time can be understood easily by looking back at Eq. dl ll i. 
Figure0b) shows this threshold temperature vs. the coupling 
constant. As expected, when the coupling constant increases 
the adiabaticity temperature-scale decreases. 



III. CONCLUSION 



^ v * M ' l + n l J (H) 

M*) = £i2(0)e ( - 2 * +K(2fí+1))í , 

where g~\\ (£12) is the diagonal (off-diagonal) element of 
g. We assume that the initial state of the spin is |^o) = 
cos f |e) + sin t\g). Unlike the case of unitary evolution of 
mixed states where the GP can be explained only by length of 
the Bloch vector (r) and the sòlid angle subtended by it on the 
Bloch sphere (Í2), i.e. $ = — arctan (r tan Ç), in the case 
of nonunitary evolution such information is not sufficient and 
the GP has a complicated dependence on the Bloch vector. 
However, dynamics of the Bloch vector can give an insight 
about how dissipative mechanisms affect the system. Figure 
[5] illustrates these effects on the Bloch vector of the density 
matrix. As is seen, in this case the length of the Bloch vector 
decreases, which characterizes a coherence loss, and, the vec- 
tor precesses toward the center of the sphere, which is related 
to dissipation. After coming back to the lab frame, to find the 
GP of the state one has to solve the eigenvalue problem of 
g(t), and plug the result into Eq. Q. The explicit form of the 



We have studied thermal effects induced by an environment 
on the geomètric phase. It has been shown that even in weak- 
coupling and adiabatic limits geomètric phase can vary with 
the temperature of the environment. One of implications of 
the results of this paper is that in geomètric phase experiments 
in systems which are in contact to a dissipative environment 
care must be taken in correct interpretation of results. The 
existence of an adiabaticity temperature-scale, similar to the 
adiabaticity time-scale, puts some constraints on the realiza- 
tion schemes of geomètric quantum information processing. 
This in turn can have some implications in the more realistic 
investigations of robustness of geomètric quantum computa- 
tion lUal. 
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